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Abstract 



We present a physical model for polarized and unpolarized structure functions 
and parton density functions (PDFs) of the proton and the neutron. It repro- 
duces the data on F|(x, Q^) for 0.00001 < j; < 1 and 2.5 < < 5000 GeV^, 
F^{x)-F^{x), F^{x)/F^{x), xg{x), d{x) - u{x), the Gottfried 

sum, the fractional momentum of charged partons and the polarized struc- 
ture functions at various Q^. We present for the first time, proton 
and neutron PDFs which do not assume charge symmetry. Contrary to the 
common practice, we explain polarized and unpolarized data with a single 
model. 
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I. INTRODUCTION 



Several new measurements of the polarized structure functions (SFs) of the proton and 
the neutron have been reported in the last few years [||]. In addition, recent experiments 
have significantly widened the kinematic range over which the unpolarized SFs and gluon 
density in the proton are known 0. The New Muon Collaboration has obtained accurate, 
final results on the ratio of the deuteron and the proton (unpolarized) SFs, which can be 
used to extract the ratio and the difference of the proton and the neutron SFs 0]. There 
are new data also on the anti-down and anti-up quarks in the proton . On the theoretical 
side, many parameterizations of unpolarized and polarized parton density functions (PDFs) 
exist 1^ . Parameter values are determined from a global fit to the relevant high-energy data. 
As new data become available, these parameterizations are generally revised. 

In this Letter, we adopt a very different approach. Our input PDFs are based on a 
different ansatz and have a few physically motivated free parameters. In principle, two of 
these parameters (a and b in (2)) could be evaluated theoretically, but in the absence of a 
full understanding of a quantum chromodynamic (QCD) bound state, namely the nucleon, 
we adopt the following pragmatic approach. We determine them by fitting data on the 
unpolarized structure function F2{x,Q'^) at only one value of Q^. Here x is the Bjorken 
variable and is the momentum scale. Thus we do not perform a global fit, with a large 
number (~ 15-20) of free parameters in the PDFs. We do not change the parameterization 
when we go from unpolarized to polarized PDFs. This is contrary to the common practice 
where they are parameterized separately. Finally, in contrast to the commonly available 
PDF sets, our PDFs for the proton and the neutron are not based on the assumption of 
charge symmetry. The issue of charge-symmetry violation (CSV) has acquired importance 
in recent months 0|. 

Once the parameters of the model are determined, it can successfully predict the rest 
of the unpolarized and polarized data, at various values of and x. We have cal- 
culated: F^{x,Q^) for 0.00001 < x < 1 and 2.5 < < 5000 GeV^, and the fol- 
lowing observables at the various values of for which data are available: -^^(a:) — 
F2{x), F2{x)/F2{x), xg{x), d{x)—u{x), (i(x)/'u(x), the Gottfried sum S'g, the longitudinal 
momentum fraction carried by quarks and antiquarks, and finally the polarized structure 
functions gi{x) and gi{x). The agreement with experimental data is nearly as good as that 
obtained by the standard PDFs 0. 

II. MODEL 

There is enough evidence to seriously consider the ansatz that the input-scale parton 
densities in the nucleon may be quasi-statistical in nature [^|[ . It is also known that statis- 
tical mechanics is applicable to isolated quantum systems with finite numbers of particles 
if the residual two-body interaction is sufficiently strong and that the interaction-driven 
statistical equilibrium that emerges in such systems can be described in terms of the usual 
statistical quantities such as temperature We do not claim that the formalism in is 
applicable in toto to the QCD bound state such as the nucleon. We do, however, believe that 
1^-0 point towards the need for having an open mind about the above ansatz. The success 
of a statistical model such as the present one, in reproducing and correlating a vast body of 
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polarized and unpolarized SF and PDF data provides a strong a posteriori justification for 
tlie above ansatz. 

The parton number density drf/dx in tlie infinite-momentum frame (IMF) and tlie den- 
sity dn I dE in tlie nucleon rest frame are related to each other by 

dn' M^x /•*^/2 dE dn 



dx 2 JxM/2 E"^ dE' 

where the superscript i refers to the IMF, M is the nucleon mass and E is the parton energy 
in the nucleon rest frame §]. This is a general relation connecting the two frames; the only 
assumption made is that of massless partons. This assumption is common in the formalism 
of deep inelastic scattering. 

There is a standard procedure in statistical mechanics to introduce the effects of the 



finite size of an enclosure, in the expression for the density of states [|10[- It allows us to 
write dn / dE for the nucleon as 

dn/dE = g f{E) {yE'^/2-n'^ + aR^E + hR), (2) 

where g is the spin-color degeneracy factor, f{E) is the usual Fermi or Bose distribution 
function f{E) = {exp[(ii^ — /i)/T] ± V is the nucleon volume and R is the radius 

of a sphere with volume V. The three terms in (2) are the volume, surface and curvature 
terms, respectively; in the thermodynamic limit only the first survives. We determined the 
parameters a and b in (2) by fitting the structure function F2{x) data at a fixed momentum 
scale Q^. Unlike a and b, the temperature (T) and chemical potential (/i) in the Fermi and 
Bose distributions are not fitted to a detailed shape of any data, but get determined due to 
number and momentum constraints on the PDFs (see below). Note also that a, b, T and yU 
are constants independent of x and Q^. 

We now depart from the procedure followed in and present a more complete model for 
polarized as well as unpolarized PDFs and SFs of protons and neutrons. If na{a)-i{i) denotes 
the number of quarks (antiquarks) of flavor a and spin parallel (antiparallel) to the nucleon 
spin, then any model of PDFs in the proton has to satisfy the following seven constraints: 

(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 

all partons 

and similarly for the neutron. The values of Au, Ad and As in (6)-(8) have been measured 
by several groups. We use Au = 0.83 ±0.03, Ad = -0.43 ±0.03, As = -0.10 ±0.03 for the 
proton and Au = -0.40±0.04, Ad = 0.86±0.04, As = -0.06±0.04 for the neutron 0. The 
summation in (9) runs over quarks, antiquarks and gluons. The numbers na(a)m) in (3)- (8) 
are obtained from (l)-(2) by integrating the appropriate dn^ /dx over x, and the momentum 
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fractions in (9) are obtained similarly by integrating the appropriate xdrt/dx over x. The 
double integrations are performed by first interchanging the order of the two integrations, 
and then doing the x integration analytically and the E integration numerically. It is nec- 
essary to distinguish between and /iaj, to ensure that the statistical model is consistent 
with the constraints (6)-(8). We note that ^la^ = —fiai and fiai = —^J^ap Hence (3)-(9) 
represent 7 coupled nonlinear equations in 7 unknowns, namely fiui, fJ^dh fJ^di, /^st' f^si 
and T. We can solve them numerically for a given choice of a and b. Once these 7 unknowns 
are determined, dn^/dx are known, and polarized and unpolarized SFs can be evaluated 
using the standard relations between SFs and PDFs. 

We have described above our model of the input- scale {Ql) PDFs. Since data are available 
at various values of Q^, it is necessary to be able to evolve these PDFs from Ql to any other 
Q^. This involves evolving singlet, nonsinglet and gluon densities in the unpolarized and 
polarized cases. We have done this by solving the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
equations in the next-to-leading order, taking = Af^, Nf = 4 and as{ml) = 0.117 |]TT|. 



III. RESULTS AND DISCUSSION 

Our fits to the structure functions -Ff (^) and -F^(x) at = 4 GeV^ are shown in Fig. 
la. We chose 4 GeV^, because the NMC data on Ff' — F2 are available only at this 
Q^. The resultant values of the two parameters, namely a and b, are given in Table I. 
We shall comment on these values later. It is easy to show analytically using (3)- (8) that 
/^mT > f^di > fJ'ui > fJ'd^ > fJ'si = —fJ'si > 0, for the proton. Figure lb shows our results 
for F2{x) — F2{x) at = 4 GeV^ in comparison with the NMC data. The discrepancy 
between the calculated results and the data, near x ~ 10~^, is of the order of 0.005 which is 
negligible compared to Ff'" in this region. The calculated Gottfried sum at = 4 GeV^, 
over the interval 0.004 < x < 0.8, is 0.215 in agreement with the experimental number 
0.221 ± 0.019(syst) ± 0.008(stat) 

Once the parameters of the model are determined, it predicts the rest of the unpolarized 
and polarized data quite successfully without any further fitting. We now present these 
results. Figure 2 shows the calculated Ff' in comparison with the data, in the range 0.00001 < 
X < 1 and 2.5 < < 500O GeV^. 

From SFs, we now turn to PDFs. The HERA and EMC data on xg{x), the Fermilab 
data on d{x) — u{x) and the CDHSW data on d{x)/u{x) are compared, in Fig. 3, with our 
results and the results based on parameterizations in In particular, the d{x)/u{x) ratio 
in the limit a; — > 1 is found to be 0.22 in good agreement with the QCD prediction 0.2 ||15|| . 
The momentum fraction carried by quarks and antiquarks, at = 15 GeV^ is 0.58 which 
is consistent with the MRST result and the experimental observation that the charged 
partons carry about half the proton momentum. 

Finally, the results on the polarized structure functions gi{x) and gi{x) in Fig. 4, show 
that the model is able to predict the shapes of the polarized data, once the parameters are 
determined as explained in Sec. II. The present model trivially satisfies the general positivity 
constraints on the polarized (Sf) and unpolarized (f) PDFs: \Sf{x,Q^)\ < f{x,Q'^). 

We sum up the differences between the present model and that presented in |g] . (a) The 
earlier version did not distinguish between /i| and /xj^. It involved solving three simultaneous 
equations, namely the two number constraints and one momentum constraint, for the three 
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unknowns, namely T, /i„ and /i^. It could not explain the polarized data. The present 
version splits up each chemical potential into two. It involves solving the seven simultaneous 
equations (3)- (9), and is able to explain and correlate the polarized and the unpolarized data 
by means of a single model, (b) We have presented here proton and neutron PDFs without 
assuming charge symmetry. Whereas the subject of charge symmetry violation in PDFs 
is being discussed extensively in the literature 0, we have presented PDFs which actually 
incorporate this feature, (c) The present model is in agreement with the observed difference 
and ratio of F|' and F^, which was not the case earlier (d) In p]], the QCD evolution 
was performed to the leading order; here it is performed up to next-to-leading order. 

Interestingly, the signs as well as magnitudes of a and h (Table I) turn out to be con- 
sistent with those in DTof, suggesting a physical basis for our parameterization. Does the 
statistical model provide a physical basis to the commonly used PDF parameterizations [^]? 
To investigate, we have compared our PDFs with those in at a common, low = 1.25 
GeV^. Our xu^ and xdy have the same shapes and similar magnitudes as those in Our 
xg is somewhat larger and xu and xd are somewhat smaller; however, their shapes are the 
same as in On this basis one would be inclined to answer the above question in the 
affirmative. Details will be published elsewhere. 

In conclusion, we have shown that the ideas from statistical mechanics work even in- 
side the nucleon. Trying to understand why they work, could deepen our understanding 
of hadron structure as well as statistical mechanics. Whereas all available PDFs assume 
charge symmetry, this paper presents, for the first time, a set of proton and neutron PDFs 
which does not make this simplifying assumption. The model is remarkably successful in 
reproducing a large body of polarized and unpolarized, PDF and SF data on the proton 
and the neutron. Thus it is potentially able to make quantitative predictions for the input 
PDFs. The scope and precision of the model can be extended systematically — e.g., by 
having a more elaborate treatment of the finite-size effects, by allowing finite mass of the 
charm quark, by considering nuclear corrections in the deuteron and higher-twist effects, 
etc. 

We are very grateful to A. Deshpande, D. Fasching and E.-M. Kabuss for many useful 
communications. 



Note added: Boros et al. proposed a large CSV of the sea quarks in the nucleon, 
as an explanation of the discrepancy between neutrino (CCFR) and muon (NMC) nucleon 
structure function data at low x. This, however, has been criticized by Bodek et al. [|] who 
showed that the above proposal is ruled out by the published CDF W charge asymmetry 
measurements. This controversy does not affect the present work because the discrepancy, 
if any, between neutrino and muon data at low x is not used as an input anywhere in the 
model. The discussion of this issue in the recent literature was used only as one of the 
motivations for this model. The origin of CSV in the present model (i.e. the origin of the 
different values of a and h for the proton and the neutron) is in the tiny differences between 
An, Ad, As values for the proton and the neutron (see the RHSs of Eqs. (6)-(8)), and also 
in the fact that we fit to and F2 separately. These two points are independent of the 
above controversy. 
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TABLE I. Values of the parameters a and b. Temperature (T) and chemical potentials (/i) are 
in MeV. T and fi are constrained by (3)-(9) once a and b are specified. 





Proton 


Neutron 


{a,b) 


(-0.376,0.504) 


(-0.300,0.504) 


T 


62 


59 




(210,86) 


(40, 94) 




(42, 106) 


(188, 76) 




(-7,7) 


(-4,4) 
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FIG. 1. (a) and F2 vs. x. The solid and dashed curves are our fits to Fj* and F2 data, 
respectively. The data are from [12,13] and the F2 data are extracted from the accurate final 
results on - at = 4 GeV^, suppHed by the NMC collaboration [3]. NMC obtained 
these results by ignoring nuclear corrections in the deuteron. Since at low x, neutron data are not 
available, the proton data were used in that region, for the purpose of the fit. (b) F2 — F2 vs. x. 
The curve represents our calculation and the data are from [3] . The small discrepancy between the 
curve and the data at low x, is discussed in the text. 
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FIG. 2. Predictions of the model, for the unpolarized structure function O'^sr a wide 

range of (GeV^) and x, compared with experimental data [2,12,13]. 
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FIG. 4. Predictions of the model, for polarized SFs. Data are from [1]. 
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